Quantum Hall states of atomic Bose gases: 
density profiles in single-layer and multi-layer geometries 
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We describe the density profiles of confined atomic Bose gases in the high-rotation limit, in single- 
layer and multi-layer geometries. We show that, in a local density approximation, the density in 
a single layer shows a landscape of quantized steps due to the formation of incompressible liquids, 
which are analogous to fractional quantum Hall liquids for a two-dimensional electron gas in a 
strong magnetic field. In a multi-layered set-up we find different phases, depending on the strength 
of the inter-layer tunneling t. We discuss the situation where a vortex lattice in the three-dimensional 
condensate (at large tunneling) undergoes quantum melting at a critical tunneling t ci . For tunneling 
well below t ci one expects weakly coupled or isolated layers, each exhibiting a landscape of quantum 
Hall liquids. After expansion, this gives a radial density distribution with characteristic features 
(cusps) that provide experimental signatures of the quantum Hall liquids. 

PACS numbers: 03.75.Lm,73.43.Cd 



I. INTRODUCTION 

Among the fascinating developments in the field of 
quantum gases is the possibility to study correlated states 
of matter in a setting that is entirely different from the 
traditional setting of electrons in a solid state environ- 
ment. A prime example are fractional quantum Hall 
states, which are expected when trapped atoms (bosons 
or fermions) are made to rotate at ultra-high angular mo- 
mentum 

HLH H II H 110 While there has been stead y 
progress in achieving high angular momentum |t| , the 
conditions for the actual realization of these states have 
not yet been met. 

The most direct experimental signature of electronic 
(fractional) quantum Hall states, the quantization of the 
Hall conductance, is not easily available for realizations 
of such states with neutral atoms. It is thus important to 
investigate the experimental signatures of atomic quan- 
tum Hall states. A number of largely complimentary 
proposals have been made 

• fractional (braid) statistics [ljj 

• loss of condensate fraction ^lj 

• detection of gapless edge excitations 

• detection of density correlations in expansion image 

m 

• characteristic density profiles 0, 0] 



In this paper, we work out the proposals put forward 
in [lj, Il5t llta ] for the detection of atomic (fractional) 
quantum Hall states via characteristic density profiles in 
a suitably engineered experimental set up. 

We first analyze the case of a single layer of rotating 
atoms, with a weak confining potential. For this situ- 
ation we find a characteristic plateau landscape for the 
density, with steps at sharply quantized values of the den- 
sity (see Figs, ^and [2J, which are direct consequences of 
the incompressibility of the quantum Hall states. 

After that we consider a more involved situation, where 
a condensate that is rotating around the z-axis is exposed 
to an optical lattice potential in the z-direction. This 
leads to a multi-layer geometry, with an inter-layer tun- 
neling whose strength is set by the strength of the optical 
lattice potential. This approach is presently pursued by 
a number of experimental groups. With an appropri- 
ate choice of parameters, one expects a regime that can 
be characterized as a stack of weakly coupled layers that 
each display a landscape of quantum Hall states. We pro- 
pose to probe this state after expansion in both the axial 
and radial directions. As the expansion of each layer in 
the axial direction is very fast, it gives rise to a density 
profile as a function of the radial distance r that is an av- 
erage over the radial density profiles of the many layers. 
Despite the fact that the density profiles in these layers 
differ from each other (owing to the differing numbers of 
particles in the layers), we show that the resulting av- 
eraged density profile has characteristic features (cusps), 
which are remnants of the steps in the individual lay- 
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ers before expansion (see Figs. 001 ■ We propose these 
cusps as possible experimental fingerprints of the atomic 
quantum Hall states. 

Our presentation is organized as follows. In section 
[n]we discuss the case of a single layer, and present the 
expected landscape of quantized plateaus. In section ITTT1 
we discuss the multi- layer set-up: section All Bl discusses 
the melting of the vortex lattice in the multi-layer case 
and section ITlI CI discusses the quantum phases expected 
after the melting. In section UlIDI we discuss the expan- 
sion after switching off the optical lattice and trapping 
potential and present the density profiles that are ex- 
pected. Throughout the paper we indicate typical values 
of the experimental parameters that are involved, and 
the conditions needed to justify the approximations that 
we make. 



II. SINGLE-LAYER GEOMETRY 

A. Experimental set-up 

We consider a single cloud of rotating bosons, with har- 
monic confinement oj±_ in the x — y (in-plane) direction 
and uju in the z (out-of-plane) direction. The character- 
istic length scales £± and £y are 

£ ul = y/h/imu^) . (1) 

We assume that the rotation frequency cu is at or near 
the critical frequency luj_, and that parameters are such 
that the system is in the two-dimensional (2D) regime. 
In this set-up the energy scale for atom-atom interactions 
is 

9qli = (2^Vi (2) 

with g = 47Th Qs , where a s is the scattering length. As- 
suming the values Hu± = 5nK, a s = 5nm, and taking 
in = 50nm for a layer defined by an optical lattice (see be- 
low), gives an energy scale g q H — 0.5nK. Since typically 
,9qH < huj^huj±, in the study of quantum Hall states, 
the single particle states can be restricted to the two- 
dimensional lowest Landau Level (LLL) Clearly, the 
observation of quantum Hall states (in single or multi- 
layer set-up) will require that all energy scales in the 
experiment (tunneling, temperature, etc.) are well below 

.9qH- 

We mostly focus on one of the following two situations, 
which we refer to as 'quadratic' and 'quartic' in-plane 
confinement. 

quadratic confinement. Here we assume a rotation 
frequency ui slightly below the critical frequency 
u>±, leaving a residual parabolic potential in the 
rotating frame of reference V^(r) = ^fc2^ 2 with 

&2 <X (u>± — U>). 



quartic confinement. This is for critical rotation u> = 
u±, but in the presence of an additional confining 
potential, which we take to be a quartic, Vi{r) = 



B. Phase separation in an external potential 

We refer to Refs. [USUI!, HUH for extensive studies 
of the quantum liquids that form at high angular momen- 
tum, after the point where the vortex lattice melts has 
been crossed. Exact diagonalization studies on a disc, 
which represent a finite-size gas in a harmonic confine- 
ment potential (for up to N = 10 particles), have revealed 
"vortex liquid" groundstates closely related to fractional 
quantum Hall states^], Q. Due to finite size effects, the 
density profiles of these states show large fluctuations. 
Features of some (but not all) of these states have been 
understood in a picture invoking non-interacting com- 
posite fermionsQ. Studies in edgeless geometries (sphere 
or torus) have established the existence of homogeneous 
quantum Hall liquids 0,3- In these geometries the liq- 
uids have a uniform density over the entire system, so 
the calculations represent regions of a trapped gas where 
the density is approximately uniform. Studies on edge- 
less geometries provide the most accurate numerical data 
concerning the bulk properties of the quantum Hall liq- 
uids that form at high rotation rate. 

In this paper we focus on inhomogeneous quantum Hall 
liquids, with density modulated by an external potential 
such as V-2(r) or Vi(r). The idea is the following. From 
the numerical studies cited earlier, we infer the existence 
of stable, incompressible quantum liquids at certain spe- 
cific filling factors Vi (the filling factor is defined below in 
terms of the ratio of the number densities of bosons and 
vortices) , and we know the energy per particle e; in these 
liquids. If we now consider a given number of particles 
in a slowly varying external potential, we can patch to- 
gether regions where the particles form specific quantum 
liquids, such as to minimize the total energy with the one 
global constraint of reproducing the total of number of 
particles in the trap. 

For inhomogeneous systems, in which the potential 
V"(r) varies slowly in space compared to £±, the density 
distribution, n(r), (the number of atoms per unit area) 
averaged on scales large compared to £± , can be obtained 
by minimizing the functional 

E = J d 2 r (e[i/(r)] + V(r)n(r) - /m(r)) , (3) 

where e[^(r)] is the interaction energy per unit area, 
which is a function of the local filling factor v(r) = 
n(r)/no [with no = l/(n£ 2 j)], and // is the chemical po- 
tential (which controls the total number of particles). 

That the interaction energy is purely local is a signif- 
icant simplification, arising from the fact that the inter- 
actions are short ranged. Minimization of the functional 
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then poses a purely local problem. At each position, the 
density n(r) is that which minimizes 

e[i/(r)] - u L n(r) (4) 

where the local chemical potential is 

fJL L ( T ) = fJL-V(T) . (5) 

For a vortex lattice (at large filling factor) , the energy 
density is 

e[u] = nobis 2 x g qH , (6) 

where b = 1.1596 is the Abrikosov parameter for a trian- 
gular lattice. In this case, the dependence of density on 
local chemical potential is simple: 

n(r) = n u L (r)/(2b) , (7) 

leading to a Thomas-Fermi profile 0] even in the LLL 

m 

For filling factors below a critical value v c (with v c on 
the order of 10 0,^1) the vortex lattice is replaced by 
groundstates that include incompressible quantum Hall 
fluids. The energy density is then a complicated function 
of v, containing cusps at the incompressible groundstates. 
The presence of these cusps gives rise to a step-like de- 
pendence of n on Taken with Eq. JSJ, this step-like 
dependence of n on /ii becomes a step-like dependence 
in space for a confined system (for which V(r) is not 
constant in space). 

We remark that non-interacting fermions in a rapid 
rotation regime, and subject to a slowly varying poten- 
tial, will display density profiles similar to the ones we 
discuss here |2J|. In that case, the result follows from 
the Landau level structure for non-interacting fermions, 
while all results discussed in this paper are due to in- 
teraction effects. The analogy between the two cases is 
strengthened by the 'composite fermion' formulation of 
some of our bosonic liquids. This analogy remains incom- 
plete, however, as some of the features that we see (such 
as the competition between composite fermion states and 
the paired Moore-Read state) cannot be captured by a 
model of non-interacting composite fermions. 

C. Example: two liquids 

As a simple example, we consider the possibility of two 
competing incompressible states: v\ = 1/2 (the Laughlin 
state, for which the interaction energy for contact inter- 
actions vanishes, e\ = 00), and vq. = 2/3 for which the 
interaction energy per particle is e 2 . These two liquids 
are the p = 1, 2 members of a series at v v = p/(p + 1), 
which can be understood as integer quantum Hall liq- 
uids of composite fermions each consisting of a boson 
with a single vortex attached 0, ■ Since the system 
is bounded, we must also consider the case of no parti- 
cles, v — 0. The energy per unit area takes the values 



e[v = 0] = e[v = 1/2] = 0, and e\v = 2/3] = (2/3)n e 2 . 
Minimizing Eq. Q with respect to these three possible 
values of the density, n — vno, we find: 

v = (fi < 0) 

v = 1/2 (0 < p, < u c ) 
v = 2/3 {fx c < n) . (8) 

with u c = 4e 2 . 

Let us take a harmonic confining potential V(r) = 
\kiT 2 , such that HL( r ) = M — \k 2 r 2 . The v — 1/2 state 
forms a disc that extends out to [il = (beyond this 
Hl < so v = 0), to a radius n = y2/i/fc2. We further 
see that there is a critical value of chemical potential, 
fi c = 4e 2 , at which the v = 2/3 state will first appear in 
the center of the trap (where fiL is maximum) . Evaluat- 
ing the total number of atoms in the v — 1/2 disc at /x c , 
we find the critical number 

A^ c = (l/2)rw? = 2e 2 /A 2 (9) 

with A 2 = ifc 2 ^. Above this critical value, the disc of 

v = 2/3 has a radius r 2 = y/2(/j, — /i c )/fc 2 . 

Thus, as the number of particles increases, at first they 
will form a disc of uniform density 1/2 no the radius of 
which increases with increasing N (/i). Once N reaches 
iV c (or n reaches n c ) a second, inner, disc of density 
2/3 no will form, and start to expand with increasing N 
(or /i). In terms of N and N c the two steps for N > N c 
are located at radii 



ri/l± = \J\{N-N c ) , r 1 /t ± = <\/~N+±N B . (10) 

The energy as a function of N takes the form 

£(A0 = A 2 (n 2 - ^6(N - N C )(N - N c ) 2 ^j (11) 

We note that the second derivative of E(N) is discontin- 
uous at N — N c . 

Note that, although we have only considered the pos- 
sibilities that v = 0,1/2,2/3 (and we have not allowed 
for intermediate values), we do not expect intermediate 
densities to appear in the density distribution. Since 
e[v] = for v < 1/2, values of v between and 1/2 
are only possible precisely at ul = and hence have no 
finite extent in a generic inhomogeneous system. For the 
range, 1/2 < v < 2/3, numerical studies indicate that 
e[v\ lies above the straight line joining the end-points 
e[v= 1/2] = and e[v = 2/3] = 2/3n e 2 . Therefore, in- 
termediate filling factors are unstable to phase separation 
into v = 1/2 and v = 2/3. 

A slightly different perspective on the formation of two 
spatially separated liquids is the following. If one imag- 
ines adding more and more particles to a v = 1/2 Laugh- 
lin droplet in a (radially symmetric) potential, one will 
eventually reach the point where the energy of an addi- 
tional particle at the edge of the Laughlin droplet exceeds 
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the energy gap Ai to create a single particle in the bulk 
Laughlin liquid at the center of the droplet. The corre- 
sponding critical value for N is N c = Ai/(2A 2 ). Using 
numerical results for Ai and e 2 (as given, for example in 
0) N c is found to be slightly higher than N c , in agree- 
ment with the fact that a v = 1/2 system with additional 
quasi-particles will lower its energy by phase separating 
into a central v = 2/3 liquid and an outer ring of v = 1/2 
liquid. 

Repeating the analysis of the co-existing v — 1/2 and 
v = 2/3 liquids in a quartic potential V±{r) 
to a critical N of 



\k±r 4 leads 



N c = 



£2 

A4 ' 



(12) 



with A4 = \k±i\, and with inner and outer edges at 



{r 2 /i A 
(n/£j 



( y/9N 2 - 8N 2 - N) , 



N -- J§N 2 - 8N 2 



(13) 




radial co-ordinate (scaled by total radius) 

FIG. 1: Density profiles in a single layer, at subcritical rota- 
tion, parametrized by density at the center of the trap. The 
results are for an energy functional obtained from exact di- 
agonalizations of a system with Nv = 6 states in toroidal 
geometry (aspect ratio a/b — \/3/2). 



The condition that V(r) be slowly varying, phrased as 

dV Jr^ l-L ^ 5qH 1 is equivalent to demanding that n/£± » 
1 for all i. For this to be possible, a total number of 
particles on the order of N = 100 or more is required. 
This puts the quantized density profiles out of the reach 
of exact diagonalization studies. 



D. More liquids 

Including in the analysis more quantum liquids, at fill- 
ing factors vi < v 2 < ... leads to a sequence of critical 
values N = Nc , marking the start of the formation of 
a central region of a quantum Hall fluid at v = v% . 

To illustrate the step-like profiles for higher filling fac- 
tors, we use the interaction energy function e[u] calcu- 
lated by exact diagonalization studies on a torus 0- 
In the graphs presented in Fig. ^ a confining potential 
V(r) = \k 2 r 2 is assumed. The energy functional used is 
for system size of Ny = 6 single-particle states on a torus 
[for which the numerical results extend up to v = 25/3]. 

For filling factors larger than v c ~ 6, the groundstate 
is a compressible vortex lattice Q • The steps seen in the 
density distribution in that regime are an artifact of the 
minimum change of filling factor, Av — 1/6, allowed by 
the numerics on a finite system with Ny =6. In this 
regime, the density profile will be smooth (we emphasize 
again that our results apply only on scales large com- 
pared to £±, so do not describe the short-range density 
modulations associated with the unit cell of a triangular 
vortex lattice). 

For filling factors less than v c ~ 6, these calculations 
show plateaus appearing at the densities of a set of in- 
compressible quantum Hall fluids, including the sequence 
of Moore- Read and Read-Rezayi states at v = k/2 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
radial co-ordinate (scaled by total radius) 

FIG. 2: Density profiles in a single layer, for cases with max- 
imum central densities of up to v — 1. The energy func- 
tional is for Nv = 12 states in toroidal geometry (aspect 
ratio a/b — l).The quantum liquid at v = 1 is the Moore- 
Read paired state and the v — 1/2 state is the Laughlin liq- 
uid. In between, a number of composite fermion states (at 
v = 2/3, 3/4) are observed. 



(black dashed lines) as well as some of the CF sequence 
v = p/(p+ 1) (black dotted lines). We note that finite 
size effects in these results for Ny = 6 can affect the 
widths and positions of the plateaus appearing in Fig. ^ 
This figure illustrates the general feature of the forma- 
tion of plateaus at the densities of incompressible states. 
One interesting detail is the competition between succes- 
sive members of the CF series at v = p/(p + 1), and the 
Moore-Read (MR) liquiddHEJ at v — 1. These states 
have very different correlations, and one does not expect 
the higher CF states to be stable against the pairing in- 
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stability that is at the basis of the formation of the MR 
state. This competition is highlighted in Fig. [3 which 
shows the density profiles for systems with maximum 
central densities of up to v = 1 calculated for a system 
with Ny — 12 vortices (which, unlike Ny — 6, also per- 
mits the appearance of v = 3/4, and for which finite size 
effects on the plateau widths are much reduced). Note 
that there is no appearance of an incompressible state at 
v = 5/6. This indicates that, at least within this small 
system calculation, the v — 5/6 state appears to be un- 
stable to phase separation into the v = 3/4 and v = 1 
states. 



III. MULTI-LAYER GEOMETRY 
A. Experimental set-up 

We now investigate the situation where an optical lat- 
tice in the z direction is imposed on a cigar shaped cloud 
of atoms rotating around the z-axis. If the optical poten- 
tial is sufficiently deep, this will define a stack of parallel 
planes with (weak) tunneling t between the planes. The 
idea is that parameters can be chosen such that, while 
the original cloud is in a mean field regime with total 
filling factor v > v c such that it would display a vor- 
tex lattice, the individual layers defined by the optical 
lattice can be in a quantum regime, so that the entire 
configuration becomes a stack of weakly coupled quan- 
tum liquids. We will assume the presence of a chemical 
potential \i(z) = /j,q — fi2{z/d) 2 (with d the distance be- 
tween the layers), which will induce a slow modulation in 
the number of particles per layer. One may view a set-up 
like this as a road towards isolating a single layer in the 
quantum Hall regime, or to perform diagnostics on the 
entire multi-layer geometry. 

Reaching the regime where quantum Hall states can 
form requires a number of conditions. First of all, to 
avoid the formation of a mean-field vortex lattice, the fill- 
ing factor in each of the layers should be sufficiently low 
(meaning high-enough total angular momentum per par- 
ticle) and the interlayer coupling t should be sufficiently 
small. The second condition is that the tunneling is well 
below the gap expected for the quantum Hall states. 

Concerning the first condition, a quick estimate is 
made as follows. The critical filling factor v c for a single 
layer has been found to be of the order v c ~ 10. Assum- 
ing N particles, Ny vortices and Nl layers, quantum 
melting will be possible if N/(Nr J Ny) < v c . Taking, for 
example, TV = 5000 particles and Ny — 100 vortices (fill- 
ing factor v = 50), slicing up the condensate into, say, 
Nl = 50 layers of 100 particles each would produce a 
filling factor per layer v' of order unity. For this choice 
of parameters, a gradual lowering of the tunneling am- 
plitude causes a quantum melting of the vortex lattice at 
a critical t = t Cl . For t < t Cl one expects a state with 
quantum liquids in each of the layers. 

The quantum Hall gap that features in the second con- 



dition is set by the interactions between the atoms (as en- 
coded in the scattering length a s ) and the length scales 
£±_ and £|| of the confinement in the x — y and z direc- 
tions. (In the multi-layer case, h refers to the thickness 
of a single layer.) A theoretical estimate^ gives 

A~0.14^.g qH ~ ^tkj x . (14) 
£ ll 

with <7 q H as in Eq. J5J. With the parameters as given 
below Eq. the quantum Hall gap is found to be in 
the order of 1 nK. 

There are several possibilities for experimental detec- 
tion of the formation of incompressible states in the 
multi-layer set-up. One possibility, which we worked out 
in [16j, is to focus on the total particle number per layer, 
i.e. on the density profile as a function of the layer co- 
ordinate p{z). In this paper we focus on a detection 
scheme where one allows a rapid expansion of the sys- 
tem and then measures the density profile as a function 
of the radial coordinate r. 



B. Quantum melting of the vortex lattice 

We now turn to the analysis of the quantum melting of 
the vortex lattice in the multi-layer geometry. For this, 
we will assume that all layers have identical axial and ra- 
dial confinement, with a chemical potential independent 
of z. If the filling factor per layer, v' = N/ (N^Ny), is 
well below the critical filling factor v c for a single layer, 
isolated single layers will be in the regime where the vor- 
tex lattice has melted. This means that for this situation, 
a gradual lowering of t (by turning up the strength of the 
optical lattice potential) , will lead to a melting of the vor- 
tex lattice at a critical tunneling t — t Cl . Equivalcntly, 
there is a critical layer filling factor v' c {t) for given tun- 
neling t. We will use two inequivalent techniques to find 
the value of v' c (t). 

Due to the short-range nature of the interaction, the 
inter-layer interactions are much smaller than both the 
intra-layer interaction and the tunneling, mean-field the- 
ory then predicts that the vortex lattices in different lay- 
ers will lie on top of each other. Note that this is very 
different from the quantum Hall multi-layer setting in 
semiconductor devices, where the long-ranged Coulomb 
interaction can dominate the tunneling. 

We model the coupled layers by including a Joscphson 
term: 

m 

where we include tunneling only between neighboring lay- 
ers, labeled by i. 

Our first approach for determining v' c {t) is to extend 
the single layer effective theory of vortex fluctuations de- 
rived by Sinova et al. [19| to the multi-layer case. We 
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FIG. 3: Critical layer filling factor v' c (t) as a function of 
ln(t/(i/(/qH)), obtained by applying a Lindemann criterion 
to the vortex fluctuations (Eq. I17H in the limit Nl — * oo, 
Ny — > oo . 

neglect the gapped sound modej^ and focus on the low- 
energy Tkachenko modes. In terms of the deviations 
u*(r) of the vortex positions from their equilibrium val- 
ues, the Josephson coupling takes the following form: 




(16) 

with an additional (constant) contribution to the chem- 
ical potential. In addition to the vortex positions, each 
layer has a phase <pi. We assume that the phases of all 
layers are identical, an approximation which is valid when 
t > t C2 . 

On distances smaller than the axial healing length (in 
units of the optical lattice period d) £ii « y/t/qquv', all 
bosons are in the same axial state. Therefore, they be- 
have as a single layer and the vortices are 'rigid' on the 
scale of In order to capture this effect, we use £|| as 
a short-distance cut-off in the axial direction. The mean 
square displacement of a vortex is given by: 

/ 2\ _ _Lv^ f rf2q 

£(M+M\e-m, ( 17 ) 

e*i y V K T V Kl J 

where Nl is the number of layers, k = 0, 2tt/Nl, ... is 
expressed in units of the inverse layer spacing, 1/d, and 
q is in units of l/(±. The constants a%, a>2 and Cqq are 
as in Ref. and 

K L (d,k) = ^ + i^i(l-cos(fc)) (18) 
q z 2 q z 

K T (q,k) = C66 q 2 + i^-i(l-COs(fc)). (19) 
2 q z 

The Lindemann criterion predicts the quantum melting 
transition to occur when the mean square displacement 




FIG. 4: Critical filling factor per layer u' c (t) as a function 
of ln(t/(u'g qH )) for iV L = 8, 50 layers and N v = 10, 100 
vortices. The solid line indicates the limit Nl — * oo, Nv —* oo. 

is|13 

(u 2 ) » a\al (20) 

where a/, w 0.1 — 0.2 is a system-dependent (a priori 
unknown) constant and ao is the lattice spacing. We 
take ah = 0.15 and have plotted the resulting critical 
layer filling factor v' c {t) in Fig. [3] 

Our second approach to estimate the critical filling fac- 
tor per layer v' c (t) is based on the depletion of the con- 
densate This criterion is, perhaps, less direct in 
determining the melting point of the vortex lattice than 
the Lindemann criterion, and may tend to overestimate 
the critical filling factor v' c (t). However, given its simplic- 
ity, it is valuable to know the results of the calculation 
based on this criterion. 

The filling factor v for bosons outside the condensate 

is: 

v-v = |^E/ B /q(W^-l), (21) 

The one-layer expressions [l^ (£°, A°) are generalized 

by &, q = £° + Kl - cos(fe)), E k ^ = ^ q - |A° I 2 . The 
k = integrand diverges when q — > 0, due to phase 
fluctuations. In a finite system we can use the system 
size as an infrared cutoff, qmin ~ 

l/^2nN v . We equate 
the quantum melting point with the filling factor where 
half the bosons are in an excited state and find the critical 
filling factors shown in Fig. 0] 

The critical filling factor for large tunneling, deter- 
mined in this way, is well approximated by v' c (t — > 00) = 
u' c (t — ► Q)/Nl, in agreement with the intuition that for 
large tunneling the system behaves as a single layer with 
total filling factor v — N L v' . As can be observed in 
Fig- EI i n a large range of tunneling amplitudes the criti- 
cal filling factor is independent of Nl and Ny- Eq. I|21|) 
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can be approximated by a converging integral, which is 
exact for Nl — > oo, Ny — > oo. For the parameters quoted 
above, Nl — 50, Ny — 100, v 1 = 1, we find that the crit- 
ical tunneling for quantum melting of the vortex lattice 
is t Cl w 0.183 qH - 



C. Quantum phases of the multi-layer systems 

The melting of the vortex lattice leaves a multi-layer 
system of highly correlated quantum liquids. The na- 
ture of the quantum liquids is not a priori clear. In the 
limit of vanishing tunneling, the isolated layers will dis- 
play single-layer behavior. The states formed close to 
the melting transition are, however, not necessarily adi- 
abatically connected to these but could be separated by 
a quantum phase transition. 

For t smaller than but rather close to t Cl there is 
the possibility for correlated multi-layer states. For the 
case of electronic quantum Hall systems, such correlated 
multi-layer states have been proposed and analyzed (an 
early reference is (24[). For the atomic systems the multi- 
layer states can be very different, as there are no inter- 
layer interactions and interlayer correlations are driven 
by tunneling alone. 

Our main interest is in the regime where the individ- 
ual layers are weakly coupled. In this regime, the atoms 
will be able to move in between layers and form an equi- 
librium density profile. In this regime we expect a local 
density approximation in both the coordinates z and r to 
be valid. In section UlI Dl below, we analyze this regime 
and work out the density profile arising after expansion. 

For t <C g q H, interlayer transport will stop and indi- 
vidual layers can be at a sharply quantized particle num- 
ber, realizing what we call a quantum Hall-Mott (qH- 
Mott) phase. To determine the boundaries of these qH- 
Mott phases, we have analyzed a simple effective Bose- 
Hubbard model. In the simplest case we take quadratic 
in-plane confinement and assume that all layers form a 
Laughlin state (with filling factor per layer v' — 1/2). 
For this situation we obtained the following expression 
for the critical tunneling t for entering the qH-Mott state 
with M particles, 



t C2 (n) = ~M[M + 



JL 

2Ac 



}[M- 



JL] 
2A 2 J 



(22) 



for a chemical potential in the interval M — | < ^ 



< 



M + |. As in the standard Bose-Hubbard model, the 
curves t c (p) define characteristic 'lobes' in the fi-t plane, 
marking the boundaries of the qH-Mott phases. We ob- 
serve that t C2 (p) has a maximum of A2/4, which is of 
order g q n/M. This scale is in agreement with the fact 
that the lowest excitations of the quantum Hall droplets 
are edge excitations with energies of order g^/M. 

In Fig.JSJwe display fJ,(z)-t diagrams for Laughlin states 
in quadratic and quartic in-plane confinement. In the 
quartic case, the maximum of t C2 grows with [i. This 




FIG. 5: Phase diagrams in the /J,(z)-t plane, for a situation 
with v = 1/2 Laughlin quantum Hall states in the layers. The 
shaded areas correspond to quantum Hall Mott states, where 
the particle number per layer is sharply quantized. In the re- 
maining part of the diagram there are interlayer fluctuations, 
giving rise to an axial superfluid. The top and bottom panels 
are for quadratic and quartic confinement, respectively. 



makes possible a situation where the central layers of 
the stack (near z = 0) display qH-Mott states, while the 
outer layers are in the fluctuating regime. 

We wish to emphasize that these quantum Hall-Mott 
phases are related to the well known Mott phase. When 
a ID lattice of atoms in the Mott phase is subjected to 
rotation (with fixed M), each site will expand to a layer 
with a well defined filling factor 1/. When v' w v c , the 
single layer critical filling factor, the system will undergo 
a quantum phase transition to the quantum Hall-Mott 
phase. This transition occurs when t < t C2 , where we 
extend the definition of t C2 to the Mott phase. 

An interesting point in the t — v' phase diagram (for 
fixed M) is (t, v') — (t*,v*), where t Cl — t C2 = t* and 
v* is the critical filling factor for t = t* . At this point, 
the Mott transition coincides with the quantum melting 
transition. 



D. Density profile after expansion 

In this section we focus on the regime where the tun- 
neling is such that we have a stack of weakly coupled 
quantum Hall layers. We can apply the ideas developed 
in section [n] where we showed how each individual layer 
is built from local pieces of quantum Hall fluid. Hav- 
ing a slowly varying chemical potential fi(z) implies that 
the density landscape varies slowly from layer to layer, 
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and that processes where atoms tunnel from the edge of 
a local quantum Hall fluid in one layer i to the corre- 
sponding edge in an adjacent layer i ± 1, are possible. 
The energy scale for such inter-layer edge excitations is 
of order gqn/Ni, with Nj the number of particles in layer 
i. We assume that the tunneling t is large on this scale, 
such that tunneling events can establish an equilibrium 
density profile n(r, z). 

The density profile can be written as 



em 



i(r, z) = n v[p - ^(z/d) 2 - V(r)} 



(23) 



with u[n] = vi for /4 < M < Pc l+1> ■ A measurement of 
particle number per layer averages this over 



2irdr n(r, z) 



(24) 



The profiles p(z) have been discussed in [T^. Here we fo- 
cus instead on the density profile arising after switching 
off the optical lattice potential in the z-direction and the 
trapping potential. In this expansion, each tightly con- 
fined layer expands enormously in the z direction, until 
all layers overlap. The expansion image is thus obtained 
by averaging over z, 



p(r) 



dz n(r, z) 



(25) 



Expansion in the radial direction, which is needed for an 
actual measurement of p{r), will not affect the density 
profile other than by an overall expanding scale|l3j. 

The integral in Eq. EH is easily worked out. If we as- 
sume quadratic in-plane confinement, V(r) — ^fer 2 , the 
profile takes the form of a sum of semi-circle arcs 



p(r) = n 



2k2 

"2 



1/2 



E 



[z/j - Z/j_i] 



with 




(26) 



(27) 



For quartic confinement the result is 



1/2 



with 



4fco - A*H 

&4 



(29) 



It is clear that, while the expansion does wash out the 
sharp steps in the density profiles, it leaves cusps in the 
averaged profiles. For in-plane confinement V(r) oc r a 
these cusps take the form of square root singularities, 
p(r) oc \/R a — r a . In Figs. Hj] and [7| we display some of 




Stf r 3 r 2 



£'1 




r 3 r 2 r L 

FIG. 6: Radial density profile p(r), for a system with quantum 
Hall liquids at v\ — 1/2, V2 = 2/3, — 3/4 and v>4 ~ 1. 
These profiles arise after expansion in the z-direction. The top 
panel is for quadratic in-plane confinement, with parameters 
chosen such that the landscape in the central layer (at z — 0) 
coincides with the top curve in Fig. The bottom panel 
assumes quartic confinement. 



the profiles p(r) . The positions of the cusps provide in- 
formation on the incompressible liquid states. Fitting the 
measured density profile to l|26|) , involving the set of fill- 
ing factors Vi, allows the critical chemical potentials, p l c , 
to be measured. This would provide direct information 
on the energetics and incompressibility of the correlated 
quantum liquids. 

We now estimate the magnitude of fluctuations that 
affect these mean field density profiles. Focusing on a qH 
fluid (vk) in one of the layers, we can evaluate the energy 
for a displacement of the outer edge, —> rk + Sr^. 
In evaluating this energy, we assume that the chemical 
potential at the center of the layer is of order g q n (which 
is the natural scale in the present regime) . With that we 
find 



SE 



(30) 



Equating this to the scale set by the tunneling gives the 
estimate 



(31) 
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FIG. 7: Radial density profile p(r), for a system with quantum 
Hall liquids at v\ = 1/2 and V2 = 2/3, for quartic in- plane 
confinement. 

We conclude that, for t of order g q n or smaller, the fluc- 



tuations induced on are of order £± or less. This means 
that, in the presence of tunneling, the steps within each 
of the layers remain well-defined on a scale £±. Accord- 
ingly, one expects that the cusps in the profile after ex- 
pansion will be recognizable if the overall extent (in the 
x-y plane) of the original cloud is large on the scale of 

ix- 

The features in p(r) can be sharpened further by elim- 
inating some of the outer layers from the stack, keeping 
only a number of layers centered around z — 0, or by 
having an axial confinement \i{z) that is steeper than 
quadratic. 
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